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ABSTRACT 

Motivated by the observed discrepancy between the strong spatial correlations of Lyman break galax- 
ies (LBGs) and their velocity dispersions, we consider a theoretical model in which these starbursting 
galaxies are associated with dark matter halos that experience appreciable infall of material. We show 
using numerical simulation that selecting halos that substantially increase in mass within a fixed time 
interval introduces a "temporal bias" which boosts their clustering above that of the underlying popu- 
lation. If time intervals consistent with the observed LBGs star formation rates of ~ 50M Q yr _1 are 
chosen, then spatial correlations are enhanced by up to a factor of two. These values roughly correspond 
to the geometrical bias of objects three times as massive. Thus, it is clear that temporal biasing must 
be taken into account when interpreting the properties of Lyman break galaxies. 
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1. INTRODUCTION 

Cosmologists love to count and correlate, and mostly 
for good reasons. If structure formation proceeds through 
gravitational instability then the number densities and 
distributions of cosmological objects are directly depen- 
dent on the underlying physical model. Thus measuring 
the number of objects as function of virial velocity and 
rcdshift provides a direct probe of the primordial power 
spectrum and the overall cosmological parameters. This 
technique has been applied most cleanly to galaxy clus- 
ters, whose number densities and evolution provide strong- 
constraints on the overall matter density (eg., Eke, Cole, 
& Frenk 1996; Bahcall, Fan, & Cen 1997; Carlberg, Yee, 
& Ellingson 1997). 

Similarly, because the density peaks in which objects 
form are more clustered than the underlying mass dis- 
tribution, the clustering of virialized structures provide a 
wealth of information. This "geometrical bias" is a system- 
atic function of the mass of these structures, an effect that 
has been well-studied analytically and numerically (Kaiser 
1984; Bardeen et al. 1986; Mo & White 1996; Porciani et 
al. 1998; Jing 1999; Scannapieco & Barkana 2002). 

Perhaps the most famous measurement of this biasing 
has been in the large sample of z ~ 3 galaxies made avail- 
able by the Lyman-break color-selection technique (Adel- 
berger et al. 1998; Steidel et al. 1998). Early papers 
pointed out that if the mass of these Lyman break galaxies 
(LBGs) could be determined, their clustering would serve 
as a sensitive test of cosmology (Mo & Fukugita 1996; 
Adelberger et al. 1998; Giavalisco et al. 1998). Later ef- 
forts inverted this approach, using a wide range of cos- 
mological constraints to compute the bias of LBGs and 
relate this to their overall mass (eg. Coles et al. 1998; Gi- 
avalisco & Dickenson 2001; Shu, Mao, & Mo 2001; Wech- 
sler et al. 2001; Porciani & Giavalisco 2002). Such com- 
parisons have shown that the clustering of LBGs brighter 
than IZab < 25.5 is roughly that expected from the geo- 
metrical bias of 10 12 M Q objects in the currently favored 
cosmological model. Furthermore, the theoretical number 
density of 10 12 M Q objects at z = 3 is consistent with the 



observed densities of LBGs, roughly 2 x I0~ 3 h 3 Mpc~ 3 , 
where h is the Hubble constant in units of 100 km s _1 
Mpc -1 . Thus counting and correlating suggest a one-to- 
one correspondence between LBGs and 10 12 M Q objects. 

Yet there are problems with this simple picture. LBGs 
are extremely luminous in their rest-frame UV, implying 
star formation rates on the order of ~ 50M Q per year (eg. 
Adelberger & Steidel 2000). Thus a one-to-one correspon- 
dence means that all 10 12 M© objects must be forming 
stars at an enormous rate at z = 3. Furthermore, the 
linewidths measured from the nebular emission of a spec- 
troscopic sample of the brightest of such galaxies gives 
projected velocity dispersions of 50 — 115 km/s (Pettini et 
al. 2001), which correspond to total masses < 10 n Mo, if 
interpreted as circular velocities. 

In this Letter we explore an alternative possibility. We 
associate LBGs with a limited subset of objects that expe- 
rience an appreciable increase in mass, which we naturally 
associate with a starburst. Through a detailed numerical 
simulation we show that such accreting groups are more 
clustered than the general population, mimicking the prop- 
erties of higher-mass halos, and modifying the mass scales 
that are most naturally associated with LBGs. 

This "temporal biasing" has never before been measured 
in simulations, and is not dependent on any merger criteria 
or properties of the accreted material. However the idea 
that LBGs correspond to merger induced starbursts has 
been proposed (eg. Kolatt et al. 1999), and several authors 
have conducted numerical studies of the bias of mergers, 
obtaining mixed results. Kauffmann & Haehnelt (2002) 
analyzed the cross-correlation between objects undergoing 
major-mergers and the general population, finding weak 
enhancement at small distances. Gottober et al. (2002) 
found that applying a merger criterion at z = 1 can effect 
the bias of objects at z = 0. Finally, Percival et al. (2003, 
hereafter P03) applied a set of merger criteria at 10 8 yr 
intervals, failing to obtain enhancement as discussed in 
detail below. 

The structure of this work is a follows: In §2 we de- 
scribe our numerical simulation and in S3 we discuss our 



1 Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5 Firenze, Italy 

2 Department of Physics and Astronomy, McMaster University, 1280 Main St. West, Hamilton, Ontario, L8S 4M1, Canada. 

1 



2 



group-finding algorithms and develop a robust definition 
of accreting groups. In §4 we present our results for the 
correlation functions of these samples, and a discussion is 
given in §5. 

2. SIMULATIONS 

Driven by measurements of the Cosmic Microwave Back- 
ground, the number abundance of galaxy clusters, and 
high redshift supernova distance estimates (eg. Spergel et 
al. 2003; Eke et al. 1996; Perlmutter et al. 1999) we focus 
our attention on a Cold Dark Matter cosmological model 
with parameters h = 0.7, Oo = 0.3, = 0.65, = 0.05, 
(Ts = 0.87, and n = 1, where £! , £1a, and Qb are the total 
matter, vacuum, and baryonic densities in units of the crit- 
ical density, erf is the variance of linear fluctuations on the 
8ft. _1 Mpc scale, and n is the "tilt" of the primordial power 
spectrum. The Bardecn et al. (1986) transfer function was 
used with an effective shape parameter of V = 0.18. 

The two competing desires of achieving high mass reso- 
lution while simulating a large sample of halos led us to use 
a box size of 73 comoving Mpc on a side, populated with 
350 3 dark matter particles. The mass of each particle was 
4.3 x 1O 8 M , which gives a nominal minimum mass reso- 
lution for our group finding of 3.4 x 10 10 M Q as we select 
only groups with 80 or more particles. The simulation was 
started at an initial redshift of z — 49, and a fixed physi- 
cal Plummer softening length of 5.7 kpc was chosen. The 
simulations used a parallel OpenMP-based version of the 
HYDRA code (Couchman et al. 1995, Thacker & Couch- 
man 2000), with 64-bit precision being used throughout. 

3. GROUP FINDING AND HALO TRACING 

Group finding is a widely studied topic in cosmology 
since it can lead to (small) systematic differences (eg 
Jenkins et al. 2001). To demonstrate the robustness of 
our results we have chosen two distinct group finding 
approaches; the friends-of-friends approach (Davis et al. 
1985, FOF) and the HOP algorithm (Eisenstein & Hut 
1998). Although it remains popular, the FOF masses esti- 
mates are known to have significant scatter due to a link- 
ing problem that can occur as small strings of particles fall 
within the linking length. 

The HOP algorithm works by using the local density 
for each particle to trace ('hop') along a path of increas- 
ing density to the nearest density maxima, at which point 
the particle is assigned to the group defined by that lo- 
cal density maximum. As this process assigns all parti- 
cles to groups, a 'regrouping' stage is needed in which a 
merger criterion for groups above a threshold density S oute r 
is applied. This criterion merges all groups for which the 
boundary density between them exceeds 5 sa ddie, and all 
groups thus identified must have one particle that exceeds 
Speak to be accepted as a group (see Eisenstein & Hut 1998 
for explicit details). 

Beginning from z — 4.89, we saved particle positions 
every 50 million years up to the final output at z = 3. For 
the final 5 outputs we found FOF groups using a linking 
parameter of b = 0.18, and HOP groups using the param- 
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eters: N dens = 48, N hop = 20, N m 
^saddle = 140, and S ou t e r = 80. Visual inspection showed 
strong similarities between the two halo populations, with 
a small amount of unavoidable noise coming from groups 



around the 80 particle resolution limit (a group found by 
FOF at this limit may not be found by HOP and vice 
versa). The group index of each particle was stored at 
each output to enable tracing between outputs. 

P03 investigated the clustering of mergers within simu- 
lations using the FOF algorithm, finding no evidence for 
bias in a series of populations defined using different se- 
lection criteria. To test their conclusions we reconstructed 
their samples within our simulation. To give a rough 
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Fig. 1. — Comparison of halo growth. The FOF algorithm 
exhibits a significant amount of scatter in mass estimates be- 
tween outputs. Only 67% of groups grow from one output to 
another, compared to 82% for HOP. 

estimation of the accuracy of the group finding methods 
in Fig. 1 we plot the mass of the most massive progenitor 
at ti(z = 3.059), versus the mass at t 2 (z = 3), such that 
At — 5 x 10 7 yr. Clearly HOP identifies groups that are 
more likely to be more massive at later outputs, while 
FOF groups show considerable scatter about the mean. 
The effect of this difference is significant. 

Using the P03 definition of 'new' groups, namely, those 
for which 50% of constituent particles were not in a pro- 
genitor of equal or higher mass at the previous time, we 
find a large difference between the total number of groups 
identified: 9367 versus 12129 for FOF and HOP respec- 
tively (from 16330 and 16831 total groups, respectively). 
However, neither of these populations has stronger clus- 
tering than the complete halo population as measured by 
their correlation functions in mass bins of 0.5 dex. Notably 
the P03 halo subsets, corresponding to 3:1 and 17:3 mass 
ratio mergers, we significantly different between the FOF 
and HOP populations. This can partially be attributed 



3 



synchronization problems, and as the time difference be- 
tween outputs is increased the identified subsets become 
more similar. Further scatter is also introduced by hard 
selection criteria, such as requiring progenitors fit a cer- 
tain mass threshold. For the 3:1 and 17:3 merger ratios 
we find only 171(226) groups between z = 3.2 and z = 3, 
respectively. The correlation functions of these halo popu- 
lations, co- added over the final five outputs, show no clear 
bias, confirming the conclusions of P03, however better 
statistics are necessary for a definitive result. 

Our definition of accreting groups is similar to the P03 
criteria, except that we select the subset of halos that grew 
by 20% from output to output, which implicitly includes 
mass accretion via smooth infall and results in 545(980) 
HOP(FOF) groups if At = 5 x 10 7 . Note that the mass 
of each halo is that at the end of each time interval, such 
that we tag all halos that experienced appreciable infall. 
The 20% value is arbitrary, but we selected it primarily 
because it appears to lie outside the central 'noise' band 
in the FOF data (see Fig 1). The number of halos corre- 
sponding to this cut is also illustrated Fig 1, as points to 
the right of the dashed lines. 

4. TEMPORAL BIAS 

In Figure 2 we show the correlation function of the 
groups selected by both the HOP and FOF algorithms and 
compare them with £(r) of the accreting groups. In the 
accreting case we co-added the correlation functions calcu- 
lated from the differences from the last four At = 5 x 10 7 
year intervals and the last two At = 10 x 10 7 year intervals. 
Radial bins of 1/80 the simulation size, corresponding to 
0.92 Mpc, were taken throughout. 

For comparison, in each panel of Fig. 2 we also show the 
correlation function of all the groups in the next largest 
mass bin. The v = 1.686D(z)- 1 a(M)- 1 values (where D 
is the linear growth factor) for each of these bins are 1.38, 
1.58, and 1.83 which result in geometrical bias factors of 
b = 1 + (z/ 2 - 1)/1.69 of 1.53, 1.89, and 2.39 respectively. 
The amplitudes of the correlation functions obtained us- 
ing all the HOP groups are in good agreement with these 
bias values, and likewise, the correlation functions of the 
full set of HOP and FOF groups agree with each other to 
within statistical uncertainties. 

The upper row of plots demonstrate a clear enhance- 
ment of the clustering of accreting groups at both the 
1O 1O 5 M0 and 1O 11O M0 mass scales, with their correlation 
functions roughly matching those of objects three times 
greater in mass (no conclusion can be drawn from the 
high mass bin as the sample is too small). This "tempo- 
ral biasing" arises from the fact that both objects accret- 
ing substructure as well as those experiencing considerable 
smooth infall tend to be found in the densest regions of 
space, which are themselves highly clustered. This con- 
clusion is supported by the fact that the average local 
overdensity of groups in the 10 10 ' 5 (10 11 )Mq mass bin is 
0.82(0.87) (measured in 4 Mpc (comoving) spheres, cor- 
responding to a mass scale of 1.2 x 10 13 M Q ), where as 
the same mass bin for the entire population exhibits an 
overdensity of 0.60(0.73). 

In the second row of Fig. 2 we take a longer interval 
of At = 10 x 10 7 yr. This has a slight dampening effect 
on temporal bias, as groups in less dense regions are able 
to sustain this level of infall. Nevertheless a detectable 



enhancement is still seen, particularly in the 1O 1O ' 5 M0 
bin, which contains the most common objects, with v the 
smallest. In the At = 20 x 10 7 yr case, however, only a 
very weak enhancement of £ (r) was measured. 

In the lower two rows of this figure, we repeat our anal- 
yses using the FOF group finder. Although this approach 
is more susceptible to statistical noise, the same trends are 
apparent as in the HOP case. If At = 5 x 10 7 yr, this tem- 
poral bias is roughly equal to the geometrical bias of the 
groups three times more massive, while if At = 10 x 10 7 
yr, £(r) is boosted to a slightly lesser degree. 

Besides the comparisons shown, we have also studied 
the effect of varying a number of parameters: S pea k, Ssaddie 
and S ou ter in the HOP method, the linking length of the 
FOF groups, and the fractional value we used to define 
"appreciable" accretion. Again, the same bias trends were 
visible, and our results remained robust over a wide range 
of reasonable choices for these quantities. 

Finally, to quantify our results, we have computed the 
effective temporal bias in each mass bin, At, and group 
finder. We define b\ as the ratio of the correlation function 
of the accreting groups to the overall correlation function, 
weighted by the number of points in each bin in the overall 
function; 6 2 = Vf „ eac " ctlng M ' Afa "-' , where the sum is car- 
ried out over all bins within r < 20 comoving Mpc. These 
values are labeled in each panel, and in the At = 20 x 10 7 
yr case, b\ = 1.1(1.0) in the 10 10 - 5 (10 n - )M© HOP bins 
and 1.1(1.3) in the respective FOF bins. 

5. DISCUSSION 

From the tests presented above, it is clear that if small 
time intervals are chosen, the clustering of accreting groups 
is robustly enhanced with respect to the underlying pop- 
ulations. To relate this to LBGs we plot the spatial corre- 
lation function of TZ ab < 25.5 LBGs, as derived by Wech- 
sler et al. (2001), in the center column of Fig. 2. Although 
there are significant uncertainties involved in computing 
this quantity, since comparisons are more naturally con- 
ducted in angular coordinates, the shaded regions provide 
a guide to the range of £(r) values consistent with obser- 
vations. In these panels, we see that if At = 5 x 10 7 yr is 
chosen, then temporal bias boosts the correlation function 
of 1O 11 M halos into reasonable agreement with observa- 
tions. 

This mass is marginally consistent with the upper mass 
bound inferred from the rotation curves of a somewhat 
bright (TZab <; 24) spectroscopic subset of LBGs (Pettini 
et al. 2001). Furthermore, only ~ 4% of all groups ex- 
hibit appreciable accretion in each At = 5 x 10 7 year time 
interval and the density of 10 11 M Q halos is ~ 2 x 10~ 2 
Mpc 3 , at z = 3 in our assumed cosmology. Thus asso- 
ciating such objects with 5 x 10 7 year starbursts results 
in a density ~ 5 x 10~ 4 Mpc 3 , comparable with that 
observed. While this time interval is small, taking the 
mean gas and dark matter ratio, and an overall star for- 
mation efficiency of 10%, results in a star formation rate 
of 0.2 x l0 n M Q n b /n o /5 x 10 7 yr = 7M Q yr" 1 . This is 
less than the observed values of ~ 50M Q yr , perhaps 
implying that an even shorter time interval is appropriate. 

While quite suggestive, these comparisons are not meant 
as a complete model, and may not prove to be the final ex- 
planation of the discrepant mass estimates of LBGs. Kine- 
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Fig. 2. — Spatial Correlation Functions. In each panel the dashed line shows the correlation function for all the groups, which the points 
connected by the dashed lines show £(r) for groups that have accreted appreciable mass in the last At years. Panels are labeled by their mass 
range and assumed At values, and in each panel, the dotted line shows the correlation function of all the groups in the next highest mass bin. 
The top two rows were generated from a set of groups selected by the HOP algorithm, while the groups in the lower two rows were selected 
using the FOF approach. The shaded region in the central panels represents the observed correlation function of IZab < 25.5 Lyman break 
galaxies as computed in Wechsler et al. 2001 by inversion of the angular correlation function. A 10% accretion threshold is applied in the 
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matic models have been explored, for example, in which 
the observed velocity dispersions of LBGs are much less 
than the circular velocities of the halos in which they are 
contained (eg. Mo, Mao, & White 1999). What is clear 
however, is that this bias can not be ignored and must 
be carefully considered when interpreting the clustering of 
these objects. While perhaps only part of the story, tem- 
poral biasing represents an important factor that must be 
taken into account when studying the properties of Lyman 
break galaxies. 
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